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Abstract: Adding a scalar singlet provides one of the simplest extensions of the Standard 
Model. In this work we briefly review the latest constraints on the mass and mixing of 
the new Higgs boson and study its production and decay at the LHC. We mainly focus 
on double Higgs production in the hh —)• bhWW —)> decay channel. This decay 

is found to be efficient in a region of masses of the heavy Higgs boson of 260 - 500 GeV, 
so it is complementary to the 46 channel, more efficient for Higgs bosons having masses 
greater than 500 GeV. We analyse this di-leptonic decay channel in detail using kinematic 
variables such as Mt 2 and the MT 2 -assisted on-shell reconstruction of invisible momenta. 
Using proper cuts, a significance of 3cr for 3000 fb ^ can be achieved at the 14 TeV LHG 
for mu = 260 ~ 400 GeV if the mixing is close to its present limit and BR(L6 —)■ hh) ~ 1. 
Smaller mixing values would require combining various decay channels in order to reach 
a similar significance. The complementarity among H —)• hh, H —?■ ZZ and H —)• WW 
channels is studied for arbitrary BR(LI —t- hh) values. 
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1 Introduction 

The discovery at the LHC by ATLAS and CMS collaborations of a scalar boson compatible 
with the Standard Model (SM) Higgs boson [1, 2] has opened a new era in particle physics. 
Since there are several Higgs production modes and five of its decay processes have been 
measured ( 77 , ZZ* ^ WW* ^ rf, and bb), it is possible to extract its couplings and compare 
with SM predictions. A useful variable to evaluate the consistency of the experimental 
data with the SM Higgs hypothesis is the so-called signal-strength modifier defined as 
fi- = O'observed/'T'SM for each channel. In the latest analyses of the full 7 and 8 TeV LHC 
data, ATLAS [3] and CMS [4] obtained following combined values 

Aatlas = i^h = 125.36 GeV), 

/icMS = 1-00 ± 0.13 {nih = 125.02 GeV) ^ ^ 

from the main Higgs production and decay modes. If we ignore the small difference in the 
ruh value used in the two fits^ and assume that there are no correlation and gaussian error, 
the combined value is given by^ 

Aatlas-i-cms = 1-10 ± 0.10. (1.2) 

^Both mh and jl should be fitted simultaneously. 

^For a recent analysis including Tevatron data and taking into account the correlation among the different 
measurements, see [5]. 
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Thus, the experimental data is certainly consistent with the SM predictions and it can be 
used to constrain new physics. 

The simplest modification to the Higgs couplings is given by a generic rescaling of 
them and this would be the case in the minimal extension of the Higgs sector, in which a 
scalar singlet that generically mixes with the ordinary SM Higgs is included. This is the 
model that we will analyze in this work. Adding a singlet to the SM scalar sector has 
implications that have been widely explored in the literature. It can help to stabilize the 
Higgs potential at high energies through their positive contributions to renormalization 
group equations that govern the Higgs quartic coupling evolution [ 6 ]. It can rescue the 
electroweak baryogenesis scenario by providing a strong enough first-order electroweak 
phase transition, as studied in refs. [7, 8 ] (see however, ref. [9]). Moreover, it can act as a 
dark matter (DM) candidate [10] or as a portal to DM [11-14], depending on its stability. 

If the new scalar is not too heavy, it can be produced at the LHC through the mixing 
with the ordinary Higgs and detected by the conventional heavy SM-like Higgs boson 
searches [15-20]. On the other hand, if the double-Higgs decay mode is open, it will 
decrease the significance of SM-like Higgs signatures. Consequently, it is important to 
explore the specific resonant double-Higgs production [21-25]. In this work we extend 
previous analysis, focusing on the particular hh —)• bbWW —)• bM'^vl~i> process, and present 
strategies to enhance the signal-background ratio by using various kinematic variables. 

The organization of our paper is as follows. In Section 2 we briefly describe the model. 
Then, we review the present constraints on the mass of the new singlet and its mixing in 
Section 3. In turn, we study the production of the heavy scalar and explore its detection 
in the next LHC run using the double-Higgs decay channel with bbi~^ui~u as a final state. 
We comment on the complementarity of this channel and the decays into two electroweak 
gauge bosons. The interplay between direct production and indirect effects, such as the 
modification of the Higgs couplings, will be considered for a luminosity of 3000 fb“^ in 
Section 4. In the next section we check the validity of our study when extending the model 
to accommodate for a DM candidate. Finally, we present the conclusions. 

2 The singlet-extended model 

One of the simplest extension of the SM Higgs sector is given by the addition of a real 
singlet field. This model has been also widely studied in refs. [21, 22, 24, 26-44]. The 
relevant Lagrangian for the scalar sector is as follows: 

C = (D^$)t(D^4.) + ^d^Sd>^S - H(4>, S), (2.1) 

with the potential [21, 29] 

F($, S) = A40 -v^y + A21 u S + A22 

+ Xo2V^S^ + Xo3vS^ + XoiS\ ( 2 . 2 ) 
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The physical doublet and singlet scalar fields can be obtained by expanding the scalar 
potential F(d>, S) around the real neutral vacuum expectation values (VEVs): 


$ = 



0 i 

+ </>/\/2 J 


S = Vs + s. 


( 2 . 3 ) 


We take v ~ 174 GeV and have chosen G(d>, S) such that vs = 0 .^ The conditions A40 > 0 , 
Ao4 > 0 , and A22 > — 2 y/X^oXoi have to be imposed in order to ensure that the potential is 
bounded from below. 


Due to the A21 term, the two scalars cj) and s mix and the mass eigenstates are given 


by 


h 


cos a sma 
— sin a cos a 


( 2 . 4 ) 


The mixing angle a and the mass eigenvalues read 


tan 2 a = 

= 


A21'\/2 
2Ao2 — A40 

(^2A4o + Ao2 ± + ( 2 A 40 — Ao2)^^ . 


( 2 . 5 ) 


The stability of the vacuum requires A02 > 0 and 4A40A02 > A21. We can use ( 2 . 5 ) to 
express (A40, A02, A21) in terms of the physical parameters a, nih, run, and v as follows: 


A40 

Ao 2 

A21 


sin^ a + mjj cos^ a 
8v‘^ 

cos^ a + mjj sin^ a 
4^2 

{mjj — m\) sin 2 a 


2 v/ 2 u 2 


( 2 . 6 ) 


The cubic and quartic interactions involving the mass eigenstates h, H can be given as 
functions of the physical parameters ( 2 . 5 ) and the three remaining couplings (A22, A03, A04). 
This is in contrast with the SM (or in the extended complex Higgs singlet model), where 
the full potential can be reconstructed from the mass (matrix) and the VEVs of the field(s). 
In what follows, we assume that h, the lighter Higgs, is the SM-like Higgs discovered at the 
LHC having nih ~ 125 GeV. Its couplings approach the SM ones in the cos a ~ 1 limit. 


3 Constraints on m/j and sin ct 

The deviation of the Higgs couplings from their SM values is constrained by the LHC Higgs 
data and by the electroweak precision observables (EWPO). We first concentrate on the 
latter. The Higgs contributes to the gauge bosons self-energies involved in the EWPO. In 

®Note that, in a generic potential, S can be shifted to fulfill this condition. 
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the singlet-extended Higgs model, the one loop self-energies will be given by the sum of two 
SM-like Higgs contributions evaluated at Higgs masses, mh and mH, rescaled by cos^ a and 
sin^ a respectively [21], This property can also be applied to non-universal diagrams {e.g., 
vertex corrections) involving the Higgses and it is transmitted to the EWPO in the limit 
where higher order, O(sin^a) terms, are neglected. Taking this into account and using 
ZFITTER [45-52], we have evaluated the predictions for the Z-peak observables [53] and 
mw, Ew [54], as a function of mn and sin^a. The list of observables used are listed in 
Table 1. The results are presented in Figure 1, where 90% and 95% C.L. allowed regions 


Observable 

Data 

Observable 

Data 

Fw 

Fz 

^had 

Re 

^FB 

80.385 ±0.015 

2.085 ±0.042 

2.4952 ± 0.0023 

41.540 ±0.037 

20.767 ±0.025 

0.0171 ±0.0010 

0.1499 ±0.0018 

sm 

Ac 

Ah 

^FB 

A0,b 

R°c 

Rb 

0.2324 ±0.0012 

0.670 ±0.027 

0.923 ± 0.020 

0.0707 ±0.0035 

0.0992 ±0.0016 

0.1721 ±0.0030 

0.21629 ± 0.00066 


Table 1. Electroweak data taken from ref. [ 55 ] used in the fit of the EWPO. 

in the niH — sin^ a plane are shown. The structure of these lines can be understood by 
noting that at mn = fnh tbe contour line is a vertical one since its value does not depend 
on the mixing angle. On the other hand, for large mu values, the mixing angle must be 
small enough to compensate the non-decoupling Higgs contributions to the EWPO. 

It is also common to use the oblique parameters (S', T, U) instead of analyzing the 
complete set of observables. We expect that in the region where mu ^ 200 GeV both 
methods should give a similar value. However, for larger mn values, the gaussian 
approximation to the that is used to fit (S', T, U) and the estimation of their errors 
starts to break down.'^ This can be explicitly checked by evaluating as a function of ruh 
using the whole Z-peak data or the oblique parameters (5, T, U). 

Let us now consider the impact of the LHC Higgs data. As already mentioned in the 
introduction, the reduction of the Higgs couplings to SM fields due to the mixing translates 
into a common reduction of the Higgs signal-strength modifier in all channels. By applying 
the narrow-width approximation, one can see that this factor is given by cos^ a. Using 
eq. (1.2), it is straightforward to derive the following bound on the mixing: sin^ a < 0.076 
(0.128) at 90% (95%) C.L. We can now combine this bound with the ones derived from 
the EWPO: the results are given in Figure 2. 

After dealing with the indirect bounds® on the mixing for a given mu value, we briefly 
comment on the direct ones, derived from heavy Higgs boson searches. Note that, as a 

^This is shown in [56], where a detailed calculation of Ar and ra^ in the singlet-extended model is 
presented. 

®There are other constraints that can be derived by imposing perturbative unitarity of scattering am¬ 
plitudes for longitudinal W bosons [38, 57[. We will ignore them since they are weaker than the other 
bounds [11]. 
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EW Precision Observables 



Figure 1. Constraints in the tuh — sin^ a derived from the full set of EWPO at the Z-peak. 

consequence of the mixing, the production and decay modes of the singlet-like Higgs H will 
be the same as those of the SM-like Higgs. However, as it has different mass, the branching 
ratios of the decay channels will be different from the SM-like Higgs. We can reinterpret 
ATLAS and CMS analyses for heavy Higgs searches to derive bounds on mu and sin^ a. 
The ATLAS collaboration presented two searches for the heavy Higgs boson. The first 
one uses H — )■ WW — )• iuiv [15] decay mode and the bound corresponds to an integrated 
luminosity of 21 fb“^ at y/~s = 8 TeV. The second one uses the H —)• ZZ decay [17]. The 
CMS collaboration has reported two analyses on heavy Higgs searches using the H —)• ZZ 
decay channels. The first one corresponds to integrated luminosity 19.6 fb“^ at ^/s = 8 
TeV and considers the i^i~qq final state [18]. The second one considers final states where 
both Z’s decay into charged leptons and corresponds to integrated luminosity 5.1 fb“^ at 
y/s = 7 TeV and 19.6 fb“^ at y/s = 8 TeV [19]. The CMS collaboration has also performed 
an analysis using the channel H — )• WW —?■ (.vtu, obtained for the conhgurations of y/s = 7 
TeV with integrated luminosity of 4.9 fb“^ and y/s = 8 TeV with 19.5 fb“^ [20]. The results 
are shown in Figure 3, where we have assumed that H has the same branching ratios as 
a SM Higgs would have for those masses. This is certainly a good approximation if the 
H ^ hh decay process is not kinematically allowed, or BR(H —)• hh) <C 1. On the other 
hand, if BR(77 —)• hh) is substantially large, these bounds have to be rescaled as indicated 
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Figure 2. Constraints in the mn — sin^ a plane derived from the full set of EWPO at the Z-peak 
combined with the LHC Higgs coupling data. We have also drawn the two benchmark points whose 
LHC implications are analyzed in detail in Section 4. 

in the figure, and eventually will become irrelevant in the BR(i? —)• hh) ~ 1 limit. In this 
case, the double-Higgs production process will be the main signature of the model at the 
LHC and deserves a detailed study. We investigate the scenario in the next section, and in 
turn, present the interplay among the different, present and future, bounds. 

4 Detecting the heavy Higgs in H ^ hh at the LHC 

The resonant double Higgs production is a distinct feature of the model we are dealing 
with. In this section we study this process in the forthcoming LHC run at 14 TeV. Since 
the Higgs production cross-section (7{H) scales as sin^ a and there is a bound on the allowed 
mixing for a given mjy, we can obtain the maximal value of (j{H) as a function mn- This 
is shown in Figure 4, where the 95% C.L. limit on sin^ a has been used. 

In order to check the signal significance at the LHC, which will be resumed with the 
upgraded center-of-mass energy along the year 2015, we perform a Monte Carlo (MC) 
study by choosing two benchmark points. For the H ^ hh decay process, the largest 
portion of signal events will consist of the four-6-jet final state as studied in refs. [58, 59]. 
However, the multi-jet signature is generically vulnerable to the huge QCD backgrounds 
and the poor reconstruction efficiency. One can attempt to increase the purity of signal 
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nifi (GeV) 


Figure 3. Constraints in the hih — sin^ a plane derived by ATLAS and CMS from SM-like heavy 
Higgs searches assuming the heavy Higgs decays as the SM one. For non-zero BR(iL —>■ hh) values, 
the vertical axis would read sin^ a/(l — BR(iL hh)). 

events by imposing a tight 6-tagging condition, but then it would result in a big sacrifice 
of the signal statistics. The ATLAS and CMS collaborations have performed searches for 
resonant double-Higgs production in the bbbb final state [60, 61]. It was shown that in order 
to be effective in this channel the mass of the new resonance should be greater than 500 
GeV to ensure two highly boosted, back-to-back 66 di-jet systems [58]. For smaller masses, 
the acceptance times the efficiency of the search decreases, thus making difficult to use this 
channel. 

The subleading decay process is bbW^W~, followed by fully-hadronic, semi-leptonic, 
and di-leptonic modes. This search channel, as it will be shown below, can be efficient in 
the range of the heavy Higgs mass 260 - 500 GeV.® Among them, the final states containing 
the lepton are more suitable for the search since the fully-hadronic states are liable to be 
in trouble due to the similar reason as in the four-6-jet signal. In leptonic signal events the 
missing energy originated from the neutrino prevents the direct reconstruction of the Higgs 
resonances. Still, provided that the light Higgs boson mass nih is accurately known, one 
can obtain the neutrino four-momenta up to a two-fold ambiguity by using on-shell mass 

®The complementary channel to the one presented here is the bbjj [62, 63]. The small branching ratio 
of the SM Higgs decaying into two photons makes this channel challenging (see, however refs. [64, 65]). 
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Figure 4. Cross-section of the double Higgs production mediated by the heavy Higgs boson 
evaluated at the maximal mixing angle. The black line shows the total cross section for this process 
while different final state cross sections are presented in colours. 

relations as well as the missing energy condition in the case of the semi-leptonic channel: 

= 0, +p^ +p‘^ + p''') = ml, (4.1) 

where is the missing transverse momenta measnred in the event, and q and q' are 
jets from the hadronically-decaying W boson. On the other hand, the on-shell relations 
are not enough to constrain the neutrino momenta in the case of the di-leptonic channel 
even though it provides cleaner signal than the semi-leptonic one. We here examine the 
discovery potential of the di-leptonic decay mode, which appears to be more challengeable 
due to the missing neutrinos although it is less vulnerable to uncertainties regarding the 
jet reconstruction, by using various kinematic variables and an approximate reconstruction 
scheme. We will see the practicability and the limitation of the search strategy in two 
different scenarios characterized by 

1. mn = 400 GeV, sin^ a = 0.06, 

2. niH = 260 GeV, sin^ a = 0.09, 

and assume BR(ff —?■ hh) = 1 for both benchmark points.^ 

^We stress that the results obtained here can be readily reinterpreted for the scenario with different 
BRjil —>■ hh) values. 
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The production cross-section is cr{gg —?■ H) x —?■ hh) = a{gg —)■ (/>) x sin^ a ~ 0.7 

(1.2) pb for niH = 400 (260) GeV in the 14-TeV LHC condition. Here, cj) is the Higgs-like 
scalar and cr{gg —)• cj)) values have been obtained from the Higgs Cross Section Working 
Group in [66] assuming that the couplings of (j) are SM-like. The search channel of interest 
is 


H ^hh^hbW^W-^h{p^)b{p^)+l+{p‘^)r{q^) + $^ (^ = e,//), (4.2) 

where the source of the missing energy is the neutrinos produced by the leptonically- 
decaying W bosons. For a numerical analysis, we have generated the MC events using 
Pythia 8 [67], interfacing with the CTIO parton distribution functions [68] for a proton- 
proton collision at ^/s = 14 TeV. The parton showering and the hadronization have been 
performed by Pythia 8 . Then, the hadron-level data has been processed with the fast 
detector-simulation program Delphes 3 [69], which reconstructs the final-state objects 
such as jets, isolated leptons, and the missing energy with the inclusion of detector res¬ 
olution effects and tagging/fake rates. The input parameters have been adjusted for the 
ATLAS detector in Delphes. FastJet 3 [70] is employed to reconstruct jets. In our sim¬ 
ulation, the anti-fcj jet algorithm [71] with distance parameter of 0.5 is chosen for the jet 
reconstruction. It is known that the tagging efficiency for the 6-jet depends on the trans¬ 
verse momentum pT and the pseudorapidity g of the jet object. Recent ATLAS and CMS 
analyses on the 6-jet identification for the experimental data indicate that the efficiency 
can be as large as ~ 60 - 80% [72]. For the sake of a simple analysis, we assume a flat 
6-tagging efficiency of 70% for pT > 30 GeV and \g\ < 2.5 and the mis-tagging efficiency is 
set to be 10% for the c-jet and 1% for the light flavor as well as the gluon jet. Isolated elec¬ 
trons (muons) are required to have pT > 13 (10) GeV within \g\ < 2.4. In order to remove 
fake leptons coming from the decays of hadrons, we discard the leptons lying within the 
angular separation ARij = ^+ Agj- < 0.4 from a jet with px > 25 GeV. Since the 
tau reconstruction efficiency is relatively poor, we reject the events containing the tau-jet 
with pt > 10 GeV. The missing transverse momentum is defined as the opposite of the 
vector sum of all the visible particles’ transverse momenta. 

Having the same final states that the signal, the di-leptonic tt process is the main 
background. The subleading backgrounds include Drell-Yan (DY), di-boson, and the SM 
Higgs processes that lead into the leptonic final states and the 6-jets. In addition, we 
consider the rare SM Higgs processes including the double-Higgs production via the gluon- 
gluon fusion (GGF), the single-Higgs production via the vector-boson fusion (VBF), and 
the Higgs boson produced in association with a weak vector boson or a top-pair, i.e., 
hWjZ and hit. The SM double-Higgs events have been generated by a modified Pythia 
6 program [73] in which the matrix element calculated with hpair [74] is implemented, 
while the other processes have been generated by Pythia 8 . We use the production cross 
section for the SM double-Higgs process obtained by hpair, which can calculate up to a 
next-to-leading order value. The tt production cross section is calculated with Top-|—h 
1.4 [75] at next-to-next-leading order, and the Higgs production cross sections except that 
of the double-Higgs process are obtained from ref. [66]. For the DY and the di-boson 
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Process 

Cross section 

H ^hh{mH = 400 GeV) 

0.66 

H ^ hh {niH = 260 GeV) 

1.18 

ti 

844.43 

GGF h 

50.35 

VBF h 

4.17 

hWjZ 

2.39 

htt 

0.61 

hh 

0.033 

DY 

91130.0 

Di-boson 

121.0 


Table 2. Production cross sections in pb for the signal and background processes at the 14 TeV 
proton-proton collision. We set rrit = 173.5 GeV and ruh = 125 GeV. 

processes, we use the leading-order cross sections calculated with Pythia 8 since most of 
them can contribute to the background by faking 6-jets and can be readily removed by 
event selection cuts, which will be discussed shortly. In Table 2, we show the cross-section 
values used in this study. 

Before going further into the analysis, we introduce one of the main kinematic variables 
and the reconstruction scheme adopted to obtain the approximate values of the invisible 
neutrino momenta. The situation with more than one invisible particle in a collider event 
is common in many extensions of the SM providing a viable DM candidate. One of the 
most studied collider variables to search for such a new physics signature is Mt 2 , which 
is a generalized transverse mass particularly known to be useful for the pair-production 
processes of new particles that eventually decay into the invisible particles [76, 77]. Suppose 
that the decay topology is like 

pp^Y + Y + U ^ V{p)x{k) + V{q)x{l) + U{u), (4.3) 

where T is a heavy unstable particle, P is a set of detectable particles such as jets or 
charged leptons, and x is the invisible particle. Here, U denotes a set of particles that 
does not participate in the decay process of Y like initial or final state radiations. For the 
new physics signature with the decay topology (4.3), the invisible momenta k and I as well 
as the particle masses my and are unknown, while only the sum of their transverse 
components can be inferred from the deficit of total transverse momentum in the collider 
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event, i.e., the missing transverse momentum. Then, Mt 2 is defined as 


Mt2 


mm 

kx + lT = prp 


max 




(4.4) 


where (f = 1, 2) are transverse masses for the decay chains, 

= my + + 2 - px • , 

=my + m\ + 2 (^-yiqxp + + m^, - qx • ly) ■ (4.5) 

Here, kx, Itj and m^ are input parameters. In practice, the transverse momenta of invisible 
particles are uniquely determined by the minimization, whereas the invisible particle mass 
rriy- is a constant that must be put by hand before the minimization. Once the correct m^^- 
value is chosen, the endpoint position of the Mt 2 distribution points to the parent particle 
mass, 


< my- (4.6) 

The situation becomes simpler when the invisible particle mass is already known as in the 
case of SM processes, where the neutrino is the only particle escaping detection and can 
be safely assumed to be massless for reconstruction purposes.® Another notable feature 
of the Mt 2 variable is that it comes in handy even when one or both parent particles are 
off-shell. This has been studied to measure the SM Higgs boson mass in the di-leptonic 
WW^*') channel [82, 83]. In the case when ruh < 2mw, at least one of the W bosons should 
be produced off-shell. Then, the maximal value of Mt 2 is not mw, but ~ m/i/2. This can 
be deduced by considering some special kinematic configurations as derived in Appendix A. 

As mentioned above, the di-leptonic system cannot be solved by on-shell mass relations 
even if Higgs boson masses are known. However, there is an approximation scheme to solve 
the unknown neutrino momenta by help of Mt 2 • When the minimization has been finalized 
to obtain the Mt 2 value, a unique solution for the transverse momentum configuration is 
picked up. One may attempt to see the correlation between these hypothetical momentum 
components and the true ones. For a subset of events whose Mt 2 values are close to M^ 2 ^, 
it can be shown that the Mt 2 solution of the transverse momenta are very close or equal to 
the true momenta. This can be justified by the fact that the Mt 2 solution is unique while 
preserving kinematic constraints,® and the endpoint of the transverse mass corresponds to 
the invariant mass of the decaying system, i.e., the parent particle mass. Then, by adopting 
the Mt 2 solution of the invisible transverse momenta in conjunction with known on-shell 
mass relations, one can calculate the longitudinal and energy components of the invisible 

®The application of Mt 2 to the SM di-leptonic tt process was firstly proposed in [78]. It is later employed 
and checked its efficiency in the real experimental analyses measuring the top qnark mass at both Tevatron 
and the LHC [79-81]. 

® The (i = 1, 2) functions are ellipses in the phase space and the Mt 2 value is determined by their 
intersecting point in the balanced configuration. This feature can be used to seek the Mt 2 value by using 
a sophisticated algorithm proposed in [84]. 
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Figure 5. (Left panel) Difference between the MAOS momentum and the true neutrino momen¬ 
tum for ruH = 400 GeV Zip'll = (p“’"“ - (Right panel) Normalized 

distributions for ran = 260 and 400 GeV and the ti backgrounds using parton-level data. 

four-momenta. This is so-called MT 2 -assisted on-shell (MAOS) approximation scheme [85]. 
One drawback of this scheme is that it cannot be applied if any of the parent particles are 
off-shell. However, it was claimed that one can circumvent the on-shell mass problem by 
substituting the transverse mass for the decay chain instead of the invariant mass into the 
on-shell mass relation [82, 83, 86]. It means that the on-shell mass relations now become 

(p + , {q + . (4.7) 

This modified scheme guarantees that there is always a real solution for the invisible mo¬ 
menta since the transverse mass is bounded from above by the invariant mass, and it 
maintains the property that the MAOS momentum is equal to the true momentum for the 
endpoint events of Mt 2 - See the left panel of Figure 5, where the efficiency of approxi¬ 
mation to the invisible momenta in the modified scheme is shown. Since the light Higgs 
boson mass here is set to be 125 GeV < 2mw, one or both W bosons produced by the 
Higgs boson are always off-shell. In this situation, the modified MAOS scheme (4.7) can be 
applied. Once the MAOS momentum has been obtained, one can construct the invariant 
mass of the total system, which corresponds to the heavy singlet-like Higgs boson mass 
given by 

= (p^ +p^+p^ + q^ + ^ ~ m^. (4.8) 

Strictly speaking, the equality holds only when = A;*''"® and The right 

panel of Figure 5, which shows distributions for the heavy Higgs signal and the SM 

double-Higgs as well as ti backgrounds using the parton-level MC event samples. One can 
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see that the peak position of the signal distribution clearly matches the mn value, while 
broad distributions are exhibited in the non-resonant background process. 

Armed with these tools, we now discuss our analysis to search for the heavy Higgs 
signal. After reconstructing the final-state objects, we select events that passed the basic 
cuts given as follows. 


• At least two isolated, opposite-sign leptons including the electron or the muon, i.e., 

e^e^, and We further require that one of them must have pT > 20 

GeV, 

• At least two 6-tagged jets with pT > 30 GeV, 

• Missing energy > 20 GeV, 

• For the opposite-sign same-flavor leptons, the event is rejected if m£+£- < 12 GeV 
to avoid the leptons produced by decays of the hadrons, and the Z-veto condition, 
which discards events containing \m£+£- — mz\ < 15 GeV, is imposed. 


We note that all the cut values have been chosen to optimize the signal significance. In 
the signal events, all the leptons are produced in the h —)• WW* decay process. In this 
case, it is known that the spin correlations of the decay mode make the charged leptons 
collinear. This feature can be used to further reduce the leptonic backgrounds. We use 
two angular cuts: the azimuthal angular difference |A(/)££| < 1.32 (1.57) and ARn = 
\J+ {A-quY < 1.34 (1.58) for the Higgs signal with niH = 400 (260) GeV. The 
upper frames in Figure 6 show clear separation between the signal and the ti background, 
particularly when rriH = 4:00 GeV. This is because the leptons can be much more boosted 
in the heavier Higgs events. The collinearity of leptons is also encoded in the other cut 
variables, the sum of the transverse momenta p^ = |p^ -|- q^| and the di-lepton invariant 
mass m£i. In the case when mu = 260 GeV, the leptons are less energetic so that the p^ 
is relatively soft. See the lower frames in Figure 6 . We require that p^ > 42 (25) GeV and 
m££ < 60 (47) GeV for the mn = 400 (260) GeV scenarios. The mu cut can also remove 
the Z — 7 - T~^T~ events in which the tau leptons decay leptonically. 

In addition to the basic selection and the leptonic cuts, one can impose cuts on the 
6 -jets. Recently, a boosted Higgs technique has been developed for processes like pp —)■ hV 
{V = W, Z) [87] or pp —)■ hh [ 88 ], followed by /i —>• 66 . In the situation where the 
Higgs boson is substantially boosted, the jets produced by the Higgs boson can often be 
considered as one fat jet, whose mass is around m/j. For very high » mh, = 

y^(A(/)bb )2 -|- [AphbY can be estimated to be 


ARhb ^ 


2m/t 

Pt 


(4.9) 


If the fat Higgs jet condition could be applied, the backgrounds, in particular, the tt events 
would be reduced very efficiently since the 6-jets in the background can have a relatively 
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Figure 6. Detector-level distributions of the kinematic variables for the two charged leptons. The 
upper frames are (Left panel) the azimuthal angular separations and (Right panel) the ARu when 
applying the azimuthal angular cut has been imposed. The lower frames are (Left panel) the sum 
of transverse momenta and (Right panel) the invariant mass mu distributions. Basic selection 
cuts are applied and all the distributions are normalized for an illustration. 


large angular separation. In the Higgs signal, p!^ can be as large as 


Pt — 


niH 


' 1 - 


4ml 


^ 156 GeV. 


m 


(4.10) 


H 


SO ARfjf, ~ 1.6 for rriH = 400 GeV in the rest frame of the heavy Higgs boson. The left 
panel in the upper frames of Figure 7 justifies this estimation. Normally, the fat Higgs jet is 
required to have AR^i, ~ 1.2 ~ 1.5 or > 200 GeV. Provided that the heavy Higgs boson 
is produced at near-threshold energy, the transverse momentum of the light Higgs has an 
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Figure 7. Detector-level distributions of the kinematic variables for the two 6-tagged jets. The 
upper frames are (Left panel) Ai?bh and (Right panel) the transverse momentum The lower 
frames are (Left panel) the di-6-jet invariant mass and (Right panel) the azimuthal angular sepa¬ 
ration between bb and £+£“ systems. Basic selection cuts are applied and all the distributions are 
normalized for an illustration. 

upper bound as given in (4.10). Therefore, we expect that the boosted Higgs technique 
will be applicable in the case of much heavier Higgs boson with mn ^ 490 GeV. 

In our benchmark points, it is inevitable to use the conventional kinematic cuts. Al¬ 
though the angular separations of the 6-jets are rather sizable, it still turns out to be smaller 
than the backgrounds in the case when mn = 400 GeV, while the cut can be applied in the 
opposite way in the case when mn = 260 GeV. This can be easily deduced from eq. (4.9), 
which predicts that the angular separation can be very large for the smaller value. On 
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the other hand, the right panel in the upper frames of Figure 7 shows that mu = 400 
GeV signal events possess much larger values of the total transverse momentum for the 
hh system. We select events with < 2.25 and > 105 GeV for uih = 400 GeV, 

while ^Rbh > 2.56 without imposing any p^ cut for mn = 260 GeV signal events. Since 
the rrih value is already known, one can further impose a cut on the di-6-jet invariant mass 
to ensure that the 6-jets are originated from the light Higgs boson. One can see that the 
invariant mass distributions have clear peaks around = 125 GeV for both benchmark 
points in the left panel in the lower frames of Figure 7. The invariant mass is required to 
have a value within a mass window of 115 (94) GeV < mbb < 146 (135) GeV for niH = 400 
(260) GeV signal. 

In the case when the heavy Higgs boson is produced near threshold, the light Higgs 
boson pair will be almost in a back-to-back configuration. Then, it is likely that the 
direction of the 66 system will be well separated from that of the system. This feature 
can be seen in the right panel in the lower frames of Figure 7, where the distributions for 
the absolute value of Acpbbji = cos“^(p^ • p^)) where px = Pt/pt, are shown. We take 
events with \A(pbb,££\ > 1-92 for rriH = 400 GeV signal. This cut is not applicable to the 
case of the mn = 260 GeV signal as the angular separation can be relatively small due to 
the small boost of each Higgs decay chain. 

We now turn to the Mt 2 cuts. For the 2b + 2i + final state, one can construct 
two kinds of Mt 2 according to the definition for the visible -|- invisible system, that is, 
either 21 -t- which contains leptons only, or 2b + 21 + which contains 6-jets as 
well as leptons when forming the visible particle system. We emphasize that Mx 2 has 
been known to be applicable to a system that can be divided into two groups of visible 
particles like processes depicted in the decay topology (4.3) with a pair production of heavy 
particles, followed by two separate decay chains. The 2i + .^x system in the signal decay 
topology (4.2) can be regarded as such a process. In what follows, the Mx 2 calculated 
for the 2£ -|- .^x system is expressed as M ^2 to distinguish it from the other kind of Mx 2 - 
As is derived in Appendix A for some kinematic configurations, the M ^2 distribution is 
bounded from above by mhj2 < mw, whereas it has a maximum at mw in the di-leptonic 
tt process since both W bosons are in on-shell. The M^2 distributions in the left panel of 
Figure 8 clearly show the endpoint structure. Another notable feature is that there are a 
number of events which have vanishing M ^2 tor both signal and background distributions. 
It corresponds to a trivial zero of M ^2 in the fully massless case, i.e., mi = my = 0 [89]. 
This happens when the missing transverse momentum lies on the smaller sector of the 
transverse plane spanned by the visible momentum vectors and q^. In this case, the 
M ^2 value is taken for a momentum partition where both transverse masses in eqs. (4.5) 
are vanishing. However, the fraction of events with the trivial zero of the Mx 2 oan be 
different depending on the preferred momentum configuration of the process. Due to the 
spin correlation and the boost, the opening angle of the charged leptons in the Higgs signal 
event is smaller than that in the di-leptonic ti events as was seen in the upper frames of 
Figure 6. It means that there are more chances to have the trivial zero of M ^2 ™ the ti 
events than in the Higgs signal. Therefore, the lower cut as well as the upper one can help 
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Figure 8. The upper frames are detector-level Mt 2 distributions for (Left panel) the 

(Right panel) 2b + 2£ + systems. The lower frame are (Left panel) and (Right panel) 

distributions for detector-level signals and backgrounds. Basic selection cuts are applied and 
all the distributions are normalized for an illustration. 

reduce the backgrounds further. This lower cut on the M ^2 increases the accuracy of 
the MAOS momenta, which will be used in the subsequent analysis. We impose the M ^2 
cut as 25 GeV < M ^^2 < 60 GeV for mu = 260 GeV signal. In the case when mn = 400 
GeV, the missing transverse momentum vector can lie inside of the opening angle of the 
di-lepton system when the light Higgs is fairly boosted. Therefore, we do not apply the 
lower cut, but only the upper cut M ^2 < 60 GeV is imposed for the mn = 400 GeV signal. 

Once M ^2 has been calculated, one can construct the invariant mass of the 2£ + 
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system by using the MAOS momentum of the invisible particle given as 


= (/ + q^ + 


(4.11) 


which equals to m/j when = /jtrue ^maos _ ^true_ j|. jg gjiown in the left panel in 

the lower frame of Figure 8. One can further employ the transverse mass of the leptonic 
system while ignoring the unknown muy value and the longitudinal momentum components 
of neutrinos, 

+ 2 ^y^IPtP + “ Pt ' 


which is bounded from above by rrih [90]. Since both distributions have distinguishable 
peak and edge structures as well as a strong correlation with m/j, we use them as cut 
variables demanding 777,™^-°® < 145 GeV and 30 GeV < < 125 GeV for mn = 400 

GeV, and 60 GeV < < 136 GeV and 58 GeV < M“ < 126 GeV for mn = 260 GeV 

signal events. We have not applied the lower cut on for mn = 400 GeV since the 

distribution is relatively distorted due to the trivial zero solutions described above. 

After counting two 5-jets as well as the charged leptons among the set of visible particle 
system, he., V = b£, one can define another kind of the Mt 2 variable, denoted as 
Recall that Mt 2 aims at the physics processes describable by the decay topology (4.3). The 
Higgs signal has a different decay topology since the invisible particle system is disjointed 
from the bb system. On the other hand, it is well known that the di-leptonic tt process is 
the one of the SM processes where the Mt 2 variable is applicable since the decay topology 
is exactly the same as (4.3), and the distribution is strictly bounded from above by 

mt- Therefore, one can still attempt to employ M^ 2 ^ to reduce the tt backgrounds if the 
edge structure of the signal distribution has a certain amount of deviation from mt- The 
M^ 2 ^ distributions for both signal and tt are shown in the right panel of Figure 8. 

Before going further, we here briefly summarize the types of the Mt 2 solutions for the 
invisible momenta. The hypothetical invisible momentum configuration that gives the Mt 2 
value can be classified in two types. One is a balanced configuration, in which 
is realized, and the other is an unbalanced one, in which ^ [77]. In each collider 

event, only one type of the momentum configuration provides the Mt 2 value, and it can 
be deduced by the invariant masses of the visible particle set in the event, my and my 
in eqs. (4.5). One can easily find that a stationary value of the transverse mass is 

attained when kx = m^px/my and lx = yirp — kx. Then, 

= my + m^, (4-13) 


which is called an unconstrained minimum of the transverse mass. Similarly, one can find 
the stationary value of = my + m^. For each stationary point, the value can 

There is an ambiguity of how to pair one 6-jet to one charged lepton because there can be two possible 
pairings in each event. Here, we define as the smaller one between two possible values of . This 

definition matches the one used to measure the top quark mass using Mt 2 in [78]. 
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(2) 

be compared to . In the case that 


Mr 


( 1 ) 




= mv + rriy > Mr 


Mt 2 is given by the unconstrained minimum of m!^\ i.e., 


( 2 ) 


lT=y*rp-kT 


(4.14) 


Mt2 = rny + rriy^- 


(4.15) 


This corresponds to the unbalanced configuration. On the other hand, if it is satisfied that 


Mr 


( 1 ) 


Mr 


( 2 ) 


\Lr^=m^p^/my 


lT=mxqT/m^ 


= mv + m^ < Mr 


( 2 ) 


= my + my < Mr 


( 1 ) 

X — 


lT=y'rp-kT 

kT=Px“*T 


(4.16) 


then Mt 2 is given by the balanced configuration in which m!^'^ = m!^\ See refs. [77, 91-93] 
for the detailed discussion of the momentum configuration types and their corresponding 
properties of Mt 2. 

In the case of M^ 2 ^ the Mt 2 value is always given by the balanced configuration since 
mi = mu = 0. On the other hand, because mu is not a constant but a variable, there exist 
sort of events in which the unbalanced configuration is selected to provide the M^ 2 ^ value. 
In the di-leptonic tt process, 


mu < 



154 GeV 


(4.17) 


when the b quark mass is neglected. Therefore, the unbalanced M^ 2 ^ has a maximum value 
smaller than mt, while the balanced Mt 2 value can be as large as m^. This means that the 
overall M^ 2 ^ distribution is bounded from above by the maximum of the balanced M^ 2 ^ 
values. For the Higgs signal, the situation is different. If one considers the case when the 
total transverse momentum of whole system is vanishing, or identically, the heavy Higgs 
has been produced at rest on the transverse plane, one can find that the balanced M^ 2 ^ 
value cannot exceed mH/‘2r by a similar consideration as done in Appendix A. However, 
by eq. (4.15), the unbalanced M^ 2 ^ has an upper bound at which can be expressed 

analytically as 


™max 

mu 


mnmw 

2mh 



~ 229 GeV 


(4.18) 


for mn = 400 GeV, while it is ~ 107 GeV for mn = 260 GeV. The maximum value in 
the above equation is achieved when one of the hypothetical neutrino momenta chosen by 
the M^ 2 ^ calculation is parallel to the momentum direction of the charged lepton sharing 
the same parent particle, while the other one is anti-parallel.^^ The mu distributions 
for various mu values and the M^ 2 ^ distributions classified by the types of the M^ 2 ^ 
solutions are shown in Figure 9, using the parton-level data for the sake of a numerical 
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Figure 9. Parton-level distributions for (Left panel) mu when mn = 260 and 400 GeV and (Right 
panel) classified by the types of the invisible momentum configuration chosen by the Mt 2 

calculation. See the text for detailed explanation. For a comparison, distributions for the di-leptonic 
tt process are shown. 

demonstration. This also means that the endpoint of as well as the mbe distributions 

for the Higgs signal events will be smaller than mt if mn < 330 GeV, and in that case, 
the upper cut instead of the lower one should be used unless the upper bound value is too 
close to mt- This observation may lead one to deduce that the efficiency of the cut 

might be the similar as that of the rubg cut. However, in our numerical study, the M^ 2 ^ cut 
turns out to perform slightly better than rribi. This might be because M^ 2 ^ incorporates 
the effect of the missing momentum and its correlation with the visible momenta. We set 
the event selection cut value as M^ 2 ^ > 165 GeV for mn = 400 GeV and M^ 2 ^ < 96 GeV 
for rriH = 260 GeV signals. 

When considering final-state particles all together, the simplest kinematic variables 
that one can construct are the invariant mass of total visible system, mhbei, and the trans¬ 
verse mass of the full system including the missing energy. Since the full visible -|- invisible 
system has a fixed invariant mass, i.e., rrifj, the invariant mass of the visible system also 
has a dependency on niH for its maximal value. One can find that 

4?77 ^ \ 

1 -^ ~ (377 GeV)^ (4.19) 

mjj J 

for mn = 400 GeV, whereas there is no definite cut-off in the ti background since is 
a variable of the event in the hadron collider. See the left panel of Figure 10. Since the 
lower bound is fixed as = mh = 125 GeV in the signal events, only upper cut on rubbii 

note that this for the Higgs signal is not a global maximum for all possible b£ pairings, but 

the maximum for a pair which leads to the smaller value of . 



- 20 







































(GeV) 


Figure 10. Normalized distributions for (Left panel) rntbee and (Right panel) distributions 

for ruH = 260 and 400 GeV and the tt backgrounds using parton-level data. 

variable can be applied. For the benchmark point with mu = 400 (260) GeV, we set the 
cut as rrihbu < 395 (200) GeV. This cut becomes important in the case of a heavy Higgs 
with lower mass value like in the case of mn = 260 GeV since it is capable of taking more 
stronger cut value. The other useful kinematic variable is the transverse mass of the full 
system defined as 

= mlbu + 2 ~ ‘ ’ (4.20) 

where + p^ + is the total visible transverse momentum. Here, the 

unknown m^i, is ignored. When all the visible particles are on the transverse plane and the 
neutrino momentum vectors are collinear so that the is vanishing, the transverse mass 
is equivalent to the invariant mass of the full system, i.e., ttih- It means that the transverse 
mass is bounded from above by mn as can be seen in the right panel of Figure 10. In real 
situation, the endpoint of the distribution is often smeared by the backgrounds and/or poor 
reconstruction efficiency of the final-state objects. Still, since the peak position is near the 
endpoint, it can provide an lower bound on mn- On the other hand, the transverse mass 
has some correlation with the MAOS invariant mass as discussed in [83]. They select the 
similar types of events in the phase space, and the efficiency is comparable to each other. 
We employ both two variables to suppress backgrounds and define the signal region. 

Combining all the cuts discussed so far, we examine their effects on the signal and 
the backgrounds by investigating how the cross sections are changing by applying event 
selection cuts. See Tables 3 and 4 for mu = 400 and 260 GeV, respectively. The unlisted 
backgrounds turned out to be almost negligible after applying the initial cuts. In summary, 
although the production cross section for mn = 400 GeV is smaller than that of mn = 260 
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Selection cuts 

H -^hh 

tt 

GGF h 

hit 

hh 

DY 

VV 

dsooo 

Basic selection 

0.54 

3560.36 

0.15 

0.072 

0.024 

272.41 

0.90 

0.48 

^Ru, Pt 

0.40 

562.02 

0.11 

0.015 

0.019 

33.56 

0.047 

0.90 

mu, Mf 2 

0.36 

314.95 

0.097 

0.009 

0.017 

11.20 

0.0 

1.1 

^maos 
nih , JVirp 

0.33 

237.96 

0.097 

0.007 

0.015 

11.20 

- 

1.2 

^Rbbi Pt 

0.23 

73.03 

0.008 

0.002 

0.012 

3.73 

- 

1.4 

rribb 

0.14 

16.24 

0.0 

~ 0.0 

0.007 

0.0 

- 

1.9 

^4>bb,U, mbbU 

0.13 

11.99 

- 

- 

0.005 

- 

- 

2.1 


0.059 

1.31 

- 

- 

0.004 

- 

- 

2.8 

Signal region 

0.048 

0.70 

- 

- 

~ 0.0 

- 

- 

3.1 


Table 3. Cut flow of signals for ruH = 400 GeV and the main backgrounds in fb. See the text 
for detailed description of the event selection cuts applied. VV denotes the di-boson processes 
(V = W, Z). (T 3000 is the signal significance calculated with a Poisson probability at 3000 fb“^ 
integrated luminosity. The signal region is defined by 345 GeV < < 425 GeV and 350 GeV < 

^maos ^ 43Q QgY 

GeV, the signal can be distinguished by several angular cut variables as well as the cut 
on M^ 2 ^. We have found a set of kinematic variables useful for the search. Eventually, 
the scenarios with a relatively lighter singlet-like Higgs boson are quite difficult to probe 
by using the kinematic event variables. In this case, one can still attempt to combine 
the search results from the other channels like bbrr and bbZZ, which have the next-to- 
subleading branching fractions, or a multivariate analysis like performed in [ 88 ]. If mn is 
much larger than 400 GeV, it is expected that the boosted Higgs technique approach is 
more promising. 

Up to now, we have assumed that BR(H —)• hh) ~ 1. This can be fulfilled in a large 
(A 22 ) Ao 3 , A 03 ) parameter-space region. We now relax this condition and suppose that 
the SM Higgs-like decays originated by the mixing are non-negligible. In this case, for a 
given ruH value, we can evaluate bounds on the mixing using the ATLAS and CMS data 
on heavy Higgs searches [16], as shown in Figure 3. The most stringent exclusion limit 
comes from the CMS search [18, 19]. This search is focused on the combination of the 
4£/2£2r final states in the H ^ ZZ channel assuming that the heavy Higgs only decays 
into SM particles, i.e., BR(H —?• hh) is vanishing. The maximal mixing angle allowed by 
this search for mn = 260 GeV is sin^ a < 0.06 (95% C.L.), while for mu = 400 GeV it is 
sin^ a < 0.11 (95% C.L.). If BR(H —)• hh) is non-vanishing, the latter constraints become 
weaker. The excluded sin^ a values for given BR(iL —)• hh) are represented in the light gray 
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Selection cuts 

H -^hh 

tt 

GGF h 

hit 

hh 

DY 

vv 

^3000 

Basic selection 

0.48 

3560.36 

0.15 

0.072 

0.024 

272.41 

0.90 

0.43 

‘^Ree, Pt 

0.28 

818.01 

0.15 

0.020 

0.022 

48.51 

0.095 

0.70 

mu, Mf2 

0.21 

206.23 

0.11 

0.006 

0.007 

0.0 

0.0 

0.80 

^maos, 

0.19 

140.69 

0.08 

0.004 

0.005 

- 

- 

0.88 


0.104 

6.65 

0.008 

~ 0.0 

~ 0.0 

- 

- 

2.21 

'fnbbu 

0.009 

3.03 

0.008 

- 

- 

- 

- 

2.82 


0.083 

2.29 

0.0 

- 

- 

- 

- 

2.99 

Signal region 

0.083 

2.19 

- 

- 

- 

- 

- 

3.06 


Table 4. Cut flow of signals for mn = 260 GeV and the main backgrounds in fb. See the text for 
detailed description of the event selection cuts applied. (T 3000 is the signal significance calculated 
with a Poisson probability at 3000 fb“^ integrated luminosity. The signal region is defined by 
180 GeV < < 265 GeV and 185 GeV < < 305 GeV. 

region in Figure 11 for both uih = 260 GeV and mu = 400 GeV. On the other hand, the 
constraints imposed by the EWPO and the LHG, shown as dark gray region in Figure 11, 
are independent of BR(Ff —)• hh). This is because they come from the modification of 
the couplings, parameterized by the mixing angle a, while the ones derived from the heavy 
Higgs searches depend directly on the value of BR(Ff —)• hh). One can also see the interplay 
between direct and indirect constraints in Figure 11. For mjy = 260 GeV, the direct search 
result on FT —)• ZZ imposes the most stringent bound, up to BR(iF —)• hh) ~ 0.4. For 
larger values of BR(Ff —)• hh), the LHG + EWPO limits are the most important since the 
direct search limit weakens. For mu = 400 GeV the direct limit is not as stringent as 
the indirect ones, which impose an upper bound of sin^ a < 0.084, independently of the 
BR(iF — 7 - hh) value. 

We can use the discovery reach of the 14 TeV LHC [94] for the Higgs boson search 
using the decay channel H —)■ ZZ —)• 41 in order to estimate the detectability of the two uih 
values as a function of the mixing and the BR(iF —)• hh).^‘^ In Figure 11, we show the 3 ct 
and 5cr significance lines for this channel for the integrated luminosity of 3000 fb“^. These 
lines show that for low values of BR(iL —)• hh) this search is able to resolve a large portion 
of the mixing angle values, leaving a small window of possible values. The sensitivity of this 
channel begin to decrease for BR(iF —?■ hh) > 0.6, just in the region where the double Higgs 
production, in particular the channel above mentioned, becomes relevant. In Figure 11, we 
have included the 3a equivalent line for the H ^ hh ^ bbWW* channel. It is important 

assume that BR(1F —>■ hh) + BRliF —>■ SM particles) = 1. 
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mn = 260 GeV 


niH = 400 GeV 


X = 3000 fb“‘ 



X = 3000 fb“‘ 



BR(H^hh) 


Figure 11. The dashed line delimits the 3cr signihcance region in the the sin^ a — BR(iX — )■ hh) 
plane for the H ^ hh ^ bhWW* 2b + 2t + 2u {^ = e, y) process for the integrated luminosity of 
3000 fb“^. The solid (dashed) black curve corresponds to the 5cr (3 ct) for the H — >• ZZ — ^ Ai/2i2T 
channels. Dark grey shaded region is the 95% C.L. CMS exclusion bounds and the light grey region 
is the one for EWPO + LHC. 

to note that both channels are complementary since they are very dependent on the value 
of BR(i7 —)• hh). As a remark the ATLAS collaboration has performed a search of heavy 
Higgses using the channel iL —>■ h/i —)• 6677 [95]. The results are not shown because the 
exclusion limit is well above the ones appeared in Figure 11. 

5 Comments on Dark Matter 

Given that the new scalar is unstable, it does not solve the dark matter problem. Never¬ 
theless, it can play a relevant role by providing a portal to DM. In this section we explore 
this possibility. The DM mass and its coupling to the new scalar will be restricted by 
requiring a DM relic density in agreement with the experimental value. We analyse the 
compatibility between this condition and the requirement of a sizable H ^ hh branching 
ratio, as assumed in the previous section. 

Let us consider an extra singlet neutral Dirac fermion transforming under a Z 2 sym¬ 
metry. There is a unique even renormalizable interaction term, so the Lagrangian gets 
enlarged by 

- mo)ilJ + XipStljil;. (5.1) 

The singlet fermion is stable due to the Z 2 parity and is then a potential, WIMP-like, DM 
candidate. 
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5.1 Relic density 


We have implemented the model in CalcHEP [96] and used the micrOMEGAs 2.4 


package [97] to evaluate the DM relic density for the two benchmark points studied in 
the previous section. The results are displayed in Figure 12, where we show the DM relic 


density as a function of the WIMP mass, m^, for different values of A^. The light red 
region corresponds to values varying from 0.001 to 1. The black solid line represents 
the relic density for A^ = 0.1. The blue band is bounded by the allowed experimental relic 
density value given by Planck [98]: 


0.1134 < < 0.1258 (95%C.L) 


(5.2) 


Note that the correct relic density can be achieved in two regions. The first one is char¬ 
acterized by a DM mass close to m/i/2, providing an enhancement of the DM annihilation 
cross section due to the resonance effect. When kinematically allowed, the Higgs decay 
into a xjj pair becomes to be dominant. As the LHC constrains the Higgs invisible width, 
that is mainly given by 



(5.3) 


this small window gets reduced (~ 1 GeV). 

There is a much wider parameter region where the enough amount of DM annihilation 
can be triggered by the heavy Higgs. Around and above the region of the heavy Higgs 
resonance, i.e., 2m^ > mu, the other annihilation channels such as V'?/: —)• hH and V’V’ 
HH are open, thus making the DM annihilation sufficient to attain the correct relic density. 
For 2m^ < rriH, the H —)> decay process will contribute to the decay width of the heavy 

Higgs boson, reducing the BR(i7 —?■ hh) ratio and then decreasing the resonant double 
Higgs production. This could affect the analysis done in the previous sector by reducing 
the statistical significance of the signal. However in the region where 2m.0 > mn the results 
would remain unaffected. For this reason we should incorporate the constraints from the 
direct detection experiments in order to know which DM regions are favoured. A similar 
study was done in [99], that agrees with our analysis. 

5.2 Direct Detection 

Direct detection experiments search for DM by means of its elastic scattering off nuclei. 
In the absence of a positive signal, present search results translate into bounds on the 
WIMP-nuclei cross section for a given WIMP mass. As the elastic scattering is produced 
at low momentum we can write the interaction as an effective operator. In our case, it is 
induced by t-channel exchange of the Higgses and is given by: 


Tleff D aq^tpi^qiqi 


(5.4) 


with [12] 



(5.5) 
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Figure 12. DM relic density as a function of the WIMP mass, for different values of A^. See 
the text for detailed description. 


The spin-independent elastic scattering cross section can be written as^^ 

SI ^ 1 r2 

where rup is the proton mass and fp is defined as 


A 

TUp 


y — 

qi=u,d,s 


+ -f^ 

+ 2jJtg 


E 

qi=c,b,t 


Qi 


Qi 


(5.6) 


(5.7) 


where the quantities fxqi represent the contributions of the light quarks to the mass of 
the proton. The full expressions for the spin-independent cross section can be found in 
refs. [12, 14]. In Figure 13 the normalized spin-independent cross section is plotted as a 
function of the DM candidate mass for the two benchmark points. This normalized cross 
section, is the product of the spin-independent cross section and the factor ^ defined 

as ^ = min{l, /0.1226}. This factor accounts for situations where V' provides only a 

fraction of the total amount of dark matter. In Figure 13 a scan over the mass and the 
Xpj parameters has been done. Only the points with a relic density equal or less than that 
from Planck are showed. The bounds imposed by LUX [101] are included as well as future 
prospect from XENON IT [102]. 

For the light DM candidate it is difficult to have the correct relic density and avoid 
the bound imposed by LUX at the same time. These conditions are compatible in a small 
region close to half of the mass of the Higgses so a resonant peak is present. However, this 
means that the decays h —)• i/u/; and H —)• i/u/; are dominant, so the model could be excluded 

^®See, for example, ref. [100] 
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Figure 13. Spin independent cross section as function of the DM mass for the two different scenarios 
with niH = 260 GeV (left) and mn = 400 GeV (right). The red line represents the bounds from 
LUX [101] while the black dashed line corresponds to the future prospects of XENONIT [102]. 

by the LHC or would spoil the results obtained in the collider analysis. Nevertheless we 
can find a region with relatively high masses of the DM candidate that fulfills both relic 
density and spin-independent cross section and is placed above the resonance produced by 
the heavy Higgs. In fact, the allowed area is induced by the opening of the ipip HH 
annihilation channel, making the cross section more effective. 

To summarize, our analysis implies that there is a region where DM requirements are 
fulfilled and is located above the heavy Higgs mass so the constraints from the LHC and 
the results obtained in the collider analysis are not affected. Furthermore we can see that 
in the next years direct detection experiments such as XENONIT are sensitive to a large 
amount of the parameter space of this model, leading to the possibility of probing it. 

6 Conclusions 

In the coming years, the LHC will further explore the properties of the Higgs boson by 
looking for possible deviation from the SM predictions [103]. In particular, after the high- 
luminosity upgrade LHC is expected to deliver 3000 fb“^ at 14 TeV [104]. This would allow 
to measure the 77 , WW^ ZZ, bb, and Higgs couplings within a 2 - 8 % error [103, 

104]. Meanwhile, the singlet-extended model is the simplest extension of the SM scalar 
sector. It predicts a universal deficit in the Higgs boson couplings to the SM fermions and 
gauge bosons caused by the mixing between the two neutral scalar states. Alternatively, 
a new contribution to the invisible Higgs width would imply the reduction of the visible 
Higgs decays, that can also be interpreted as a generic Higgs coupling deficit. The direct 
production and detection of the new Higgs would certainly elucidate this point. Since 
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the relevant cross section depends on the mass and the mixing, we have first reviewed 
the present experimental bounds on these two parameters. Concerning the constraints by 
EWPO, we have improved previous analysis by using the full set of electroweak observables 
instead of the oblique parameters (S', T) since the last ones only provide an accurate 
descriptions of the heavy Higgs effects in the mn ~ region. 

In order to illustrate the detection of the direct heavy Higgs production, we have chosen 
two benchmark points compatible with present bounds, in particular, the LHC Higgs data 
and the EWPO. We have studied the resonant SM Higgs boson pair production in the 
hh —)• bb WW —7- hbi^vi~v decay channel. The main background to the signal is the 
di-leptonic tt process. Besides some basic selection cuts, we have applied Mt 2 cuts for the 
2i + or 2b+ 21 + systems to optimise the signal significance. Using the di-leptonic 
channel alone, a significance 3a for 3000 fb ^ can be achieved at the 14 TeV LHC for 
mjj = 400 GeV if the mixing is close to its present limit and BR(iL —)■ hh) ~ 1. A 
lower branching ratio or a smaller mixing angle would require combining various hh decay 
channels. The complementarity between H ^ hh and H ^ ZZ channels is studied for 
arbitrary BR(LI —)• hh) values. 

We have also checked that it is possible to extend the model by including a DM 
candidate. The next generation of direct detection experiments will be capable of probing 
a large amount of the parameter space of the model. 

Note added: After completion of this work, some similar results have been presented 
in [105]. 
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Appendix 

A Higgs AfT 2 in di-leptonic WW process 

We here consider the maximum of Mt 2 in the di-leptonic decay mode of the h —)• WW^*'> 
process, 

pp ^ h + j ^ WW^*^ + J —t (-{p)y{k) + i{q)u{l) + j{u), (A.l) 
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where j denotes the initial state radiation, typically jets in the final state. The transverse 
components of the total momentum should be conserved and therefore 

Px + kx + (jx + lx T ux = 0. (-^-2) 

Since the visible particle in each decay chain, i.e., the charged leptons are massless, the 
Mx 2 value is always achieved in a balanced configuration in which m!^'^ = or 

(IPtI + |kx |j -|px + kx I — (|qx| + dx 1) “IQt + Ix I > 

X + lx = —pt — qx — ux, (A-dj 

where and 1 ™®'°® stand for the Mx 2 solution for the invisible transverse momenta. 

The above equations are satished when = —qx and 1™®'°® = —px for vanishing ux. 

On the other hand, if ux is sizable, the solution can be redefined as 

k”-- = -q^ - ^ - 5 ut 

iros = -PT-^ + ^ux. (A.4) 

where (5ux is a function parameterizing the transverse boost effect of the solution by ux 
while preserving conditions (A.3) and (5ux(ux = 0) = This solution is generically 
different from the true invisible momenta, i.e., kip®"”® ^ kx and 1 ™®"°® ^ lx, however by 
construction, the sum of each component must be equal, 

j^maos ^ imaos ^kx +lx. 

The maximum value of Mx 2 for given visible momenta and the sum of invisible mo¬ 
menta can be deduced from a kinematic property as given below, 

Mx < ruh, (A. 6 ) 

where Mx and mh are the transverse and the invariant masses for two charged leptons and 
two neutrinos dehned as 

M| = + 2 y^|px+ k “'^°®|2 y|qx + kp'"°T+ (a^x^)^ 

- 2 (px + kpi"°®) • (qx + 1 ^“), 

'ml ={p + k + q + lf, 

respectively. The transverse mass can be fu 
M^ = -\C\‘^ + B‘^ 

where 

A = IpxI + |kfp^°® 

C = PX + k“^°®. 

See [89] for the dedicated discussion of the Mt 2 solution in the case of fully massless visible and 
invisible particles. 


(A.7) 


'ther simplified as 
- |D|2-h2^5-2C.D, 


i? = |qx| + |li 

D = qx + f 


maos 

•T I 


maos 

T 


(A. 8 ) 

(A.9) 
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One can see that the equations in (A.3) is satished for the kinematic configuration where 


B = A, D = eC (e^ = 1). 

(A.IO) 

Here, e = — 1 corresponds to a kinematic configuration with vanishing ux, while it is 
non-vanishing and dux = px — qx in the case of e = 1. Then, it is true that 

M| = 4A2 - 4|C|2 - (e - l)(e + 3)|C|2 > 4^^ - 4|C|2 = 4 

(A.ll) 

for both cases. Since Mx 2 = = M^\ the relation (A. 6 ) implies that 


Mx 2 < 

(A.12) 

where the equality holds when eqs. (A. 10) are satished. Up to now, we have 
assumption whether the parent W boson is on-shell. In the case that both 
on-shell, Mx 2 is bounded from above by mw- Collectively, one can express 
of Mx 2 as 

not made any 
W bosons are 

the maximum 

Mt 2 < min jmiy, . 

(A.13) 
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